Classification of the electronic correlation strength in the Fe-pnictides: 
The case of the parent compound BaFe2As2 
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Electronic correlations in the Fe-pnictide BaFe2As2 are explored within LDA+DMFT, the com- 
bination of density functional theory with dynamical mean-field theory. While the correlated band 
structure is substantially renormalized there is only little transfer of spectral weight. The com- 
puted k-integrated and k-resolved spectral functions are in good agreement with photoemission 
spectroscopy (PES) and angular resolved PES experiments. Making use of a general classifica- 
tion scheme for the strength of electronic correlations we conclude that BaFe2As2 is a moderately 
correlated system. 
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The recent discovery of an entirely new class of high- 
Tc superconductors based on quasi-twodimensional FeAs- 
layers Q rather than on CuO-layers has stimulated in- 
tense experimental and theoretical activity. In striking 
similarity with the high-Tc cuprates, the undoped pnic- 
tides are not superconducting under ambient pressure, 
but exhibit an antiferromagnetic commensurate spin- 
density wave below 150 K [2]. When electrons (or holes) 
are added to the system via doping antiferromagnetism 
is suppressed and superconductivity appears. Since it 
is generally accepted that the strong Coulomb correla- 
tions between the Cu 3c? electrons are responsible for the 
anomalous properties of cuprates, it is tempting to sug- 
gest that the same is true for the Fe 3(i electrons in pnic- 
tides. 

Investigations of the Coulomb correlation strength 
in the Fe-pnictides have arrived at different conclu- 
sions. Haule, Shim and Kotliar Q employed the 
LDA-I-DMFT approach [4], the combination of density 
functional theory (DFT) in the local density approxima- 
tion (LDA) with the many-body dynamical mean-field 
theory (DMFT) 0, to study LaOi_:rF:j;FeAs. Using a 
value of the local Coulomb interaction U = A eV ob- 
tained by the random phase approximation (RPA) Q 
for metallic iron Q (which is quite similar to that found 
by constrained LDA, U = 3-4 eV [1, H, they con- 
cluded that LaOi_a:Fa;FeAs is a strongly correlated, bad 
metal which is close to a Mott metal-insulator transition. 
A similar conclusion was reached by Craco et al. for 
SmOi_a:Fa:FeAs. By contrast, using constrained RPA 
to compute U in LaFeAsO Nakamura, Arita and Imada 
[l^ obtained smaller values of U , in the range U — 2.2- 
3.3 eV. Therefore they concluded that in LaFeAsO and 
LaFePO electronic correlations are moderately strong. 
This was affirmed by Anisimov et al. ^ who performed 



LDA-t-DMFT investigations for LaFeAsO, with U taken 
from constrained LDA , [l3 | , and compared with pho- 
toemission spectroscopy (PES) data. Yet another con- 
clusion was reached by Yang et al. who investigate 
the Fe-pnictides by XAS and RIXS and compared with 
results from cluster diagonalizations, multiplet and DFT 
calculations. They estimated the Coulomb interaction as 
[/ < 2 eV and inferred that the Fe-pnictides are weakly 
correlated systems. 

These different assessments raise a general question 
about the classification of the strength of electronic cor- 
relations in a particular material. An overall estimate is 
provided by the ratio of the local Coulomb interaction U 
and the band width W . For U /W < 1 the system is con- 
sidered weakly correlated, and the results of DFT approx- 
imations are sufficient to explain its electronic and mag- 
netic properties. By contrast, if U is comparable with W 
or larger, the system is moderately or even strongly cor- 
related, and the Coulomb interactions must be treated 
explicitly in electronic structure calculations. However, 
the single parameter U /W is only a rough measure of 
the electronic correlations; furthermore it is not directly 
accessible by experiment. Here the effective mass renor- 
malization m* /m of quasiparticles is a more useful pa- 
rameter since it permits experimental tests. Nevertheless 
a single number cannot reliably assess the correlation 
strength in a real material whose properties are deter- 
mined by various different aspects such as the electronic 
band structure, degeneracies, filling, and the Coulomb in- 
teraction parameters beyond U. A much more expressive 
quantity is the k-resolved single-particle spectral function 
yl(k, oj), and its k-integrated variant A{lu), both of which 
can be measured by angular resolved photoemission spec- 
troscopy (ARPES) and PES, respectively. They allow 
for an unambiguous classification of the strength of elec- 
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tronic correlations. Namely, if A(k, ui) is well described 
by DFT approximations it is justified to call the system 
weakly correlated. If the band structure is substantially 
renormalized by correlation effects as expressed by an en- 
hancement of the quasiparticle effective mass, but A(uj) 
does not yet show incoherent Hubbard bands, the system 
is moderately correlated. Only if the correlation induced 
transfer of spectral weight is so strong that pronounced 
Hubbard bands and a distinct quasiparticle peak at the 
Fermi energy appear, can the system be called strongly 
correlated (a good example is the 3d^ system SrVOa (16]). 
The spectral functions of correlated electron materials 
can be computed within the LDA+DMFT framework ^ , 
which provides access to material specific single-particle 
spectra and higher correlation functions. 

In this Letter we employ the LDA-I-DMFT scheme Q 
to explore the importance of electronic correlation effects 
in BaFe2As2, the undoped parent compound of one of 
the main pnictide high- Tc superconductor families. For 
this purpose we calculate the spectral functions A{u!) and 
A(k, to) and compare our results with the available ex- 
perimental data for BaFe2As2. In particular, we apply 
the scheme described above to classify the correlation 
strength in this important class of Fe-pnictides. 

The LDA-I-DMFT approach proceeds as follows: First, 
an effective Hamiltonian is constructed using con- 
verged LDA results for the system under investigation, 
then the many-body Hamiltonian is set up, and finally 
the corresponding self-consistent DMFT equations are 
solved. By projecting onto Wannier functions [l^ we 
obtain an effective 16-band Hamiltonian which incorpo- 
rates five Fe d orbitals and three As p orbitals for each 
of the two Fe and As ions per unit cell. In the present 
study we construct Wannier states for a single energy 
window including both p and d bands. Consequently, the 
eigenvalues of the effective Hamiltonian exactly cor- 
respond to the 16 Fe and As bands from LDA. Thereby 
hybridization effects between As p and Fe d electrons are 
explicitly taken into account. The LDA calculations were 
performed with the experimentally determined crystal 
structure [18[ using the Elk full-potential linearized aug- 
mented plane- wave (FP-LAPW) code [l^. Parameters 
controlling the LAPW basis were kept to their default 
values. The calculated LDA band structure €lda{^) is 
found to be in good agreement with that of Fink et al. 
0. 

To account for the Coulomb interaction already 
present in LDA the dd-diagonal elements of the effective 
Hamiltonian are renormalized by a double counting 
correction f3| Edc = C/(n-dmft - \), 
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Here ndmft is the total, self-consistently calculated num- 
ber of d electrons per Fe site obtained within the 
LDA-j-DMFT. This form of Edc yields rehable results for 



transition metal compounds, including the superconduc- 
tor LaFeAsO Q . The p — d Hamiltonian to be solved by 
DMFT then has the form 



h'^^ d] Vu-wrr + hf^ pI dh 



(2) 



Here dkao- and pk7<T are Fourier transforms of diaa and 
Pi-ym which annihilate the d or p electron with orbital 
and spin indices aa or in the iih unit cell, and nf^^ 
is the corresponding occupation number operator. 

The DMFT self-consistency equations were solved it- 
eratively for imaginary Matsubara frequencies. The aux- 
iliary impurity problem was solved by the Hirsch-Fye 
quantum Monte Carlo (QMC) method [U. The ele- 
ments of matrix are parameterized by U and J 
according to procedure described in In the fol- 
lowing we use the interaction parameters J7=3.1 eV and 
J=0.81 eV obtained in [l3|. Calculations were per- 
formed in the paramagnetic state at the inverse temper- 
ature (3 = l/r=20 eV~^ . The imaginary time interval 
< r < /3 was divided into 150 slices, and 6x10^ QMC 
sweeps were used in a self-consistency loop within the 
LDA-fDMFT scheme. The local self-energy Ts{lu), which 
is formally an 16x16 matrix with the only nonzero ele- 
ments on the diagonal of the dd block, was calculated for 
real energies uj by analytic continuation from imagi nary 
Matsubara frequencies using the Pade approximant [23| . 

A quantitative measure of the electron correlation 
strength is also provided by the quasiparticle renormal- 
ization factor Z = {1 — ^\uj=q)'^ which leads to an 
effective mass enhancement m*/m = Z~^. In general, 
the self-energy is a matrix, leading to different effective 
masses for different bands. The masses can be obtained 
from the self-energy on the real axis, The calcu- 

lated m* /m values for every d-orbital are presented in Ta- 
ble 1. The (1^2 _y2 orbital has the smallest effective mass 
renormalization (m*/m = 1.83). The other d-orbitals 
have a slightly larger value {m* /m — 2.05 - 2.07). The 
overall effective mass enhancement m* /m Ri 2 agrees well 
with the ARPES data for BaFe2As2 [mHI. It is also in 
good agreement with experimental results for the pnic- 
tide material LaFePO, where de Haas- van Alphen exper- 
iments [26] found a mass enhancement between 1.7 and 
2.1, and ARPES [S^] measured a band renormalization 
by a factor of 2.2 compared with the LDA. 

The orbitally resolved Fe Sd and total As Ap spectral 
functions computed within the LDA and LDA-I-DMFT, 
respectively, are compared in Fig. [T] The tetragonal 
crystal field splitting of the Fe-3d states in this mate- 
rial is rather weak (Ac/=0.25 eV). Within the LDA all 
five Fe d orbitals form a common band in the energy 
range (—2, -1-2) eV relative to the Fermi level (band width 
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TABLE L Effective mass renormalization m* /m of quasipar- 
ticles in BaFe2As2 for different orbitals of the d shell. 
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FIG. 1: (Color online) Orbitally resolved Fe 3d and total As p 
spectral functions of BaFe2As2 obtained within LDA+DMFT 
(solid lines) are compared with LDA results (shaded areas). 



W ^ A eV) . There is a significant hybridization of the Fe 
t2g orbitals with the As p orbitals, leading to spectral 
weight in the energy interval (—3, —2) eV where the As 
p band is located. However, this is not a correlation ef- 
fect. The correlations only lead to some broadening but 
not to substantial spectral weight transfer, i.e., there are 
no Hubbard bands found in the entire spectrum. Indeed, 
the overall shape of the LDA spectral function is hardly 
changed by the correlations. 

In Fig[2]the total Fe 3d and As 4p spectral functions 
are compared with photoemission spectroscopy data of 
de Jong et al. [28] . The computed spectra reproduce the 
experimental features very well, including the relative po- 
sitions of the d and p bands and the strong d contribution 
at the top of the valence band. 

We now calculate the k-resolved spectral function 

A(k, u) = -Im-Tr[[uj + ^)/ - k - (3) 

TT 

Here is the 16x16 Hamiltonian matrix on a mesh of 



•a 




-5 -4 -3 

Energy (eV) 

FIG. 2: (Color online) Normalized total Fe d and As p spectral 
functions of BaFe2As2 calculated within LDA-I-DMFT (upper 
curves) are compared with photoemission spectroscopy data 
of de Jong et al. ;28l] (lower curves). 
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FIG. 3: (Color online) The k-resolved total spectral func- 
tion A(k, Lo) of BaFe2As2 near the F and X points in the 
Brillouin zone is depicted as a contour plot. Upper panel: 
LDA-fDMFT spectral function including the renormalized 
band structure (circles) obtained by plotting the peak po- 
sitions of the spectral function A(k,t<;). Lower panel: The 
corresponding experimental ARPES intensity map of Liu et 
al. '2%\. 
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k-points and /i is the self-consistently determined chem- 
ical potential. In Fig. [3] we compare our results with 
ARPES data of Liu et al. [2^. Both theory and ex- 
periment show dispersive bands crossing the Fermi level 
near the F and X points. In addition, two weakly pro- 
nounced relatively flat bands can be seen in the region 
from -0.2 to -0.3 eV and from -0.5 to -0.6 eV near the 
F point. The calculated shape and size of the electron 
and hole pockets centered at the X and F points, re- 
spectively, are in good agreement with the ARPES data. 
Also shown in Fig. [3] (upper panel) is the correlated 
band structure eoMFrO^)- Near the Fermi energy, i.e., 
in the energy range from -0.25 eV to zero where quasi- 
particles are well defined (as expressed by a linear behav- 
ior of ReS(a;)), this dispersion is very well represented 
by the scaling relation eoMFTO^) = eLDA(k)/(m*/m), 
with m* /m as the computed mass enhancement. This 
shows that the band structure is substantially renormal- 
ized. Nevertheless, since there is no substantial spectral 
weight transfer in the spectral function, electronic corre- 
lations in BaFe2As2 are only moderately strong. 

In conclusion, by employing the many-body dynam- 
ical mean-field theory in combination with LDA elec- 
tronic band structure calculations (LDA+DMFT ap- 
proach) we calculated the single-particle k-integrated 
and k- resolved spectral functions for BaFe2As2, a 
parent compound of the superconducting iron pnic- 
tide family. The correlated band structure obtained 
within LDA+DMFT is well described by the relation 
eDA/_FT(k) — eLDyi(k)/(m*/m), with an effective mass 
enhancement computed as m* /m « 2. At the same time 
correlations do not lead to substantial transfer of spec- 
tral weight in the spectrum, i.e., Hubbard bands do not 
form. Our results are in good agreement with PES ex- 
periments of de Jong et al. [28| and ARPES data of 
Liu et al. 29]. We also discussed a general classification 
scheme for the electronic correlation strength based on 
the renormalization of the band structure and the trans- 
fer of spectral weight. Our results obtained for BaFe2As2 
and their comparison with experimental data show that 
this material, and also other Fe-pnictides, should be clas- 
sified as moderately eorrelated. 
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